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We show that a kink and a topologicaUy trivial soliton in the Gross-Neveu model form, in the 
large- limit, a marginally stable static configuration, which is bound at threshold. The energy of 
the resulting composite system does not depend on the separation of its solitonic constituents, which 
serves as a modulus governing the profile of the compound soliton. Thus, in the large-A*' limit, a 
kink and a non-topological soliton exert no force on each other. 
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The problem of finding the particle spectrum (e.g., 
bound states) of quantum field theory is a major ob- 
jective of nonperturbative studies thereof. This issue 
may be addressed quantitatively in model field theories 
in 1 -I- 1 space-time dimensions such as the Gross-Neveu 
(GN) model 0], in the large N hmit. The GN model (and 
other similar models) are particularly appealing, since 
they exhibit, among other things, asymptotic freedom 
and dynamical mass generation, like more realistic four 
dimensional models. 

One version of writing the action of the 1 4- 1 dimen- 
sional GN model is 
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where the ■(/'a (a — 1, ■ • • , N) are N flavors of massless 
Dirac fermions, with Yukawa coupling to the scalar aux- 
iliary field a{x). This action is evidently symmetric under 
the simultaneous transformations a — > ~a and tjj — > 75^/', 
which generate the so-called discrete (or K2) chiral sym- 
metry of the GN model. The GN action has also flavor 
symmetry 0{2N), which can be seen by breaking the N 
Dirac spinors into 2N Majorana spinors. Related to this 
is the fact that the model is also invariant under charge- 
conjugation 0. Thus, focusing on bound states of the 
Dirac equation 



(2) 



associated with the GN action, if ipb{x,t) — e *'^''*Uf,(a;) 
(with < < m^) is a bound state solution of (jSJ, 



so is its charge conjugate spinor ilj^{x,t) 



u^ix), 



SO that bound state frequencies come in pairs: ±cjf,. If, 
however, Q has a bound state at uji, = 0, it is of course 
unpaired, i.e., self-charge-conjugate. 

Performing functional integration over the grassman- 
nian variables in the GN action leads to the partition 
function Z = J Da exp{iS'e//[(T]} where the bare effec- 
tive action is 



Seff[cr] = - 
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and the trace is taken over both functional and Dirac 
indices. 

The theory Q has been studied in the limit N 
00 with Ng"^ held flxed 1]. In this limit the partition 
function Z is governed by saddle points of Q and the 
small fluctuations around them. (In this Letter, as in 
0, we will consider only the leading term in the 1/A^ 
expansion, and thus will not compute the effect of the 
fluctuations around the saddle points.) The most general 
saddle point condition reads 
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In particular, the non-perturbative vacuum of the GN 
model is governed by the simplest saddle point of the 
path integral associated with it, where the composite 
scalar operator ipij) develops a space-time independent 
expectation value, signaling the dynamical breakdown 
of the discrete chiral symmetry by the non-perturbative 
vacuum. Thus, the fermions acquire mass m dynamically. 

Associated with this breakdown of the discrete symme- 
try is a topological solitoujthe so called Callan-Coleman- 
Gross-Zee (CCGZ) kink |l H 0, a{x) = TOtanh(m2;), 
with mass M^mk = is the dynamically generated 

fermion mass) . It is topology which insures the stability 
of these kinks: they are the lightest topologically non- 
trivial solitons in the GN model. The Dirac equation Q 
in the kink background has a single self-charge-conjugate 
bound state at = 0, which can populate at most N 
valence fermions. Thus, it gives rise to a multiplet of 2^ 
degenerate states, with mass equal to Mkink, which can 
be identified as the (reducible) spinor representation of 
Oi2N) 0,0. The expectation value of the fermion num- 
ber operator Np in a state in which the kink traps n va- 
lence fermions is n — where the subtracted piece is the 
so-called "fractional part" of the fermion number |^ |^ , 
due to fluctuations of the fermion fleld in the topologi- 
cally nontrivial background. Thus, in the kink multiplet 
-f < A^F < f . 

The GN model bears also non-topological solitons. 
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which were discovered by Dashen, Hasslacher and Neveu 
(DHN) [3 (after the work of Henceforth, we shaU 
refer to them as "DHN sohtons" . These non-topological 
sohtons are stabihzed dynamicallyjby trapping fermions 
and releasing binding energy. In la, DHN used inverse 
scattering analysis [3 to find static soliton solutions to 
the large-A^ saddle point equations of the GN model. 
(DHN also found in oscillatory, time dependent so- 
lutions of the saddle point equations, which we will not 
disucss in this Letter.) The remarkable discovery DHN 
made was that all the physically admissible static, space- 
dependent solutions of I0J, i.e., the static bag configu- 
rations in the GN model (the CCGZ kink being a non- 
trivial example of which) were reflectionless. That is, the 
static cr(a;)'s that solve the saddle point equations of the 
GN model (subjected to the obvious boundary condition 
ct(±oo) = ±m) are such that the reflection coefficient of 
the Dirac equation ((SJ associated with the GN action, 
vanishes identically. 

The Dirac equation (O in a DHN soliton background 
has a pair of charge conjugate bound states at ztujt,. The 
0{2N) flavor symmetry mixes particles and antiparticles. 
At the level of the Dirac equation |5J this means that we 
have to consider the pair ztw;, of bound state eigenfre- 
quencies together in the following way: Due to Pauli's 
principle, we can populate each of the bound states 
zbcjf, with up to N (non- interacting) fermions. Then all 
the multi-particle states in which the negative frequency 
state is populated hy N ~ h fermions (i.e., has h holes) 
and the positive frequency state contains p fermions, with 
h+p — n fixed are degenerate in energy, and thus form a 
C'^\ = {2Ny./n\{2N-ny. dimensional irreducible 0(2N) 
multiplet, namely, an antisymmetric tensor of rank n 0- 
Superficially, < n < 2N, in accordance with Pauli's 
principle. However, for dynamical reasons, as explained 
in Section 3 of 8] , only solitons with < n < are real- 
ized. The expectation value of the fermion number oper- 
ator Np in a state in which the DHN soliton traps p parti- 
cles and h holes is simply Np = p—h = 2p—n, i.e., purely 
the naive valence contribution. (There is no fracional 
contribution due to the trivial topology.) Thus, in the 
DHN sohton multiplet -n < Np < n. DHN found that 
in this case lui, — m cos ■ The mass of such a soliton 



is Af„ = sin (1^) and its profile is a{x) ^ a{oo) 
K tanh [k{x - 
where a{oo) 



«(x-.o) + ilog(^) 



±TO, K 



xa)] — Ktanh 
— cr(— oo) = 
and xq is a translational collective mode. 
Note that both M„ and Mkink are of order iV ~ p-, 
as typical of soliton masses in weakly interacting QFT. 
The binding energy Bn = nm — M„ of the DHN soli- 



ton (as well as the binding energy per-fermion, Bn/n) 
increase with the number n of trapped valence fermions 
(this is the so-called "mattress effect" in the physics of 
fermion bags). Thus, a DHN soliton is stable against 
decaying into a bunch of non-interacting fundamental 
fermions. It is also stable against decaying into lighter 
DHN bags, because any such presumed process can be 
shown to violate either energy or fermion number con- 
servation. Thus, DHN solitons are stable. Note that 
Mn=N — = 2Mkink, and more over, that a,s n N 
(i.e., K — > m and u^b 0), the profile (t(x) tends to 
cr(cx))-|-TO tanh (rnx) — m tanh ['m{x + R)] , i? — > oo . Thus, 
the configuration at n = A is that of infinitely separated 
CCGZ kink and anti-kink bound at threshold. This sin- 
gular behavior occurs because a DHN soliton, being a 
topologically trivial configuration, cannot support a nor- 
malizable zero mode, as explained below. 

The soliton solutions (both topological and non- 
topological) in the GN model serve as concrete calculable 
examples of fermion-bag |^ ^3 formation. Furthermore, 
these "multi-quark" bound states of the GN model are 
analogous to baryons in QCD in the limit of large num- 
ber of colors Since the work of DHN, these fermion 
bags were discussed in the literature several other times, 
using alternative methods [l^ . For a recent review on 
these and related matters, see [l^ . 

In this Letter, we show that the spectrum of the 
GN model contains, in the large-A^ limit, a composite, 
marginally stable topological soliton, which may be in- 
terpreted as a kink and a non-topological soliton bound 
at threshold. Furthermore, the energy of this system does 
not depend on the intcrsoliton distance, which thus serves 
as a modulus controlling the shape of the corresponding 
static solution of the large- A^ saddle point equation. 

It is a general feature of the Dirac equation in 
the background of a static topologically non-trivial a{x) 
configurations, that the spectrum contains an unpaired 
bound state at lu — 0. For example, for kink bound- 
ary conditions {a{oo) = — cr(— oo) = m), such a normal- 
izable zero-mode is given by uq cxp — cr{y) dy, where 

ij^UQ = -Uq. 

To make our point, we have to find a topologically 
non-trivial static solution of l@J, which in a certain limit 
appears as well separated kink and a DHN soliton. Thus, 
we must find a reflectionless ct{x) configuration, such that 
(O has a bound state at ui = 0, as required by topol- 
ogy, and a pair of charge conjugate bound states at some 
zbcjfc 0, with Ub considered a free parameter. To this 
end, we apply the machinery of (see also jlj]) and find 
the most general such reflectionless background as 
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2(m + k) 



1 + ^e2m(x-xo) 
(m — k)^ 



\m — K J \m — K J \m — K J 



K tanh [k{x — j/o + R)] 



sinh [m{x — xq) + k{x — yo) + 2kR] + sinh [m{x — xq) — k{x — yo)] 
g-Kfl cosh [771(0; — xq) + k{x — yo) + 2kR] + e"^ cosh [77i(x — Xq) — k{x — yo)] 
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where k = J vn? — w? , nR = i loe; 221+^ and where 
a;o and j/q are arbitrary real parameters. The latter 
two quantities arise in the inverse scattering formal- 
ism as arbitrary parameters, independent of the param- 
eter LOb, which determine the coefficients in front of 
the asymptotic exponentially decaying normalized bound 
state wave functions. Thus, the normalized bound state 
at w = behaves asymptotically as •\/2m exp —m{x—xo), 
and the normalized bound states at ±Wb behave asymp- 
totically as -\/2Kexp— k(x — yo), as x 00. Here, 
of course, xo and yo are translational collective coordi- 
nates of the soliton 17(0;). Note that 17(3:) in JSJ) satisfies 
kink boundary conditions: a{oo) = m and cr(— 00) = 
777-|-2k — 2(«;-|-77i) = —777. Similarly, ~(7(x) is the desired 
extremal configuration with boundary conditions of an 
anti-kink. 

The effective action Seff in (|3J), evaluated at the back- 
ground Q, is an ordinary function Sef f{i^b,yo — xo) of 
the parameters which determine the shape of (O, where 
we have invoked translational invariance of Q. (With 
no loss of generality, we can always set one of these col- 
lective coordinates, say xo, to zero.) Minus the value 
of Sef f{u!b,yo — Xo) per unit time is the rest energy, or 
mass M{uJb,yo — xo), of the static inhomogeneous con- 
densate ©■ We still have to extremize M{ujb,yo — xo) 
with respect to ujb- As in the case of the DHN soliton, 
the extremal value is determined by the number n of va- 
lence particles and antiparticles which are trapped in the 
bound states iwf,. Following the same technique used by 
DHN in 0] to calculate the mass of the DHN soliton, we 
find that the extremal value is again uJb = m cos (|^), 
and that the mass of ©, which we will refer to as the 
"heavier topological sofiton" (HTS), is 
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Thus, MHTS,n coincides with the sum of masses of 
a CCGZ kink and a DHN soliton trapping n valence 
fermions. Clearly, the 0{2N) quantum numbers of the 
HTS are those of the direct product of the 2^ dimen- 
sional spinorial representation and the antisymmetric 



tensor representation of rank 77. More details of the con- 
struction of lO and the associated extremum condition 
on ojb will be given elsewhere 

Note that MHTS,n is independent of the remaining col- 
lective coordinate yo- By varying yo (while keeping Ub 
fixed at its extremal value), we can modify the shape 
of cr{x) in (0) without aflFecting the mass of the soliton. 
The translational collective coordinate yo is thus a flat 
direction of the energy functional, or a modulus. In the 
following we will show that the modulus yo is essentially 
the separation between a CCGZ kink and a DHN soliton 
which we interpret as the loosely bound constituents of 
the HTS. Thus, the fact that dMHTS,n/dyo = 0, means 
that these solitons exert no force on each other, whatever 
their separation is. 

This is a somewhat surprising result, since one would 
normally expect soliton-soliton interactions to be of the 
order p- ~ iV in a weakly interacting field theory, which 
is consistent, of course, with what one should expect from 
general counting rules. Indeed, drawing further the 
analogy between the solitons discussed in this Letter and 
baryons in QCD with large N color, the HTS would corre- 
spond to a dibaryon. From the general counting rules 
[ill, the baryon-baryon interaction is expected to be of 
order N . Yet, due to dynamical reasons which ellude us 
at this point, the solitonic constituents of the HTS avoid 
these general considerations and do not exert force on 
each other. 

It is straightforward to obtain the asymptotic behavior 
of ||SJ) for large \yo\ (and xo =0). In the limit yo —00, 
the shape of Q is that of a CCGZ kink, with a little 
"DHN bump" , centered on its left wing at Xbump = J/o + 
(3/4k) log ^^^^^^ ~ yo, of width 1/2k and maximum 
value of 777 — 2ajf,. Its shape is given approximately by 
"■(xbump -I- 2;) — —m + 2k^ jim + ujb cosh2Kz). In the 
other limit, yo +00, jS)) has the shape of a CCGZ 
kink, with a little "DHN dip" , centered on its right wing 

at xdip = yo + (1/4k) log yo, of width 1/2k 

and minimum value of 777 — 2^(777 + k — LOb)/{m + k + 
UJb). Its shape is given approximately by (T(xdip + z) ~ 
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m — 2K{m + k) / {m + K + ujb exp 2kz) + 2KUJb/{c^b + {m + 
k) exp2Kz). For yo in the range such that K\yo\ ~ 1, the 
"DHN disturbance" and the kink partly overlap. These 
statements are demonstrated in the figures, for the case 
m — 2k — 1. 



Y 




FIG. 1: The soliton © for m = 2fi: = 1; 2nyo = -20 




FIG. 2; The soliton © for m = 2k = 1; 2Kyo = -1 

The limit n ^ is of some interest. Strictly speak- 
ing, there is no HTS with n ^ N, since a.t n = N the 
pair of bound states at ±uJb of the Dirac equation Q 
would coincide with the bound state at w = 0, which al- 
ready exists due to non-trivial topology. Clearly, such a 
degeneracy cannot occur in the spectrum in one spatial 
dimension. (For more details see Sections 3.1.2 and B.3.2 
in 0.) However, it is possible to study HTS's with n ar- 
bitrarily close to N . In this case, k ^ m, with i? ^ oo in 
(O. Thus, for |a;|, \xo\, \yo\ << R, ® tends in this limit 



FIG. 3; The soliton (gj for m = 2k = 1; 2K,yo = 4-20 

to 

I _ g-2m(a;-j/o) _ ^-2m{x-xo) 

a(x) = ni- — ^ -—. r (7) 

(where for clarity of presentation we have reinstated the 
parameter xq into the expression for (j(x)). In the asymp- 
totic region 1 << to|xo — j/ol (<< "mFC) Q simplifies 
further, and appears as a kink mtanh[m(a; — Xmax)], lo- 
cated at a;niax = max{a;o, ?/o}- This clearly has mass 
Mkink = but according to ©, MHTS,n~N should 

tend to = 'iMkink- The extra mass 2Mkink corre- 

sponds, of course, to the kink-anti-kink pair which re- 
ceded to spatial infinity. 

Finally, we must settle the important issue of stabil- 
ity of the HTS. Due to conservation of the topological 
charge q = (ct(oo) — (t(— oo))/2m, the final static con- 
figuration will obey the same kink boundary conditions 
as (jSJ. The HTS is too light to decay into a configura- 
tion of CCGZ kink-kink-antikink (see ©). Thus, it can 
only decay into a CCGZ kink plus a bunch of lighter 
DHN bags and/or free fermions. The binding energy 
Bn = nm + Mkink — MHTS,n of the HTS coincides with 
that of a DHN bag with the same quantum number n. 
Thus, similarly to the latter, the HTS is stable against 
evaporation into a CCGZ kink and a cloud of non inter- 
acting fermions. Decay of the HTS (with 0{2N) quan- 
tum number n ) into a CCGZ kink and a DHN bag (of 
quantum number n') is almost totally forbidden: Energy 
conservation obviously requires < n' < n. Fermion 
number conservation, on the other hand, requires (see 
for more details) n < n' < N , i.e., the complimentary set 
of the possible range (0, N) for n'. The two conservation 
laws are compatible only at rt' = n. Thus, the mass of 
the decay products equals the mass of the parent HTS, 
and the allowed channel in this case has no phase space. 
Finally, one can show, using elementary considerations as 
above, that energy conservation and fermion number con- 
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servation strictly forbid decay of the HTS into a lighter 
HTS and a DHN soliton, or into a CCGZ kink plus any 
number of DHN solitons, or into any final state contain- 
ing the time dependent solitons discovered by DHN in Q . 
Since the only allowed channel has no phase space, the 
HTS is marginally stable. This must be a manifestation 
of the fact that the translational collective coordinate yo 
is a flat direction of the energy functional. 
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